Quantum reactive scattering of ultracold NH(X ) radicals in a magnetic trap 
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We investigate the ultracold reaction dynamics of magnetically trapped NH(X '^E") radicals using rigorous 
quantum scattering calculations involving three coupled potential energy surfaces. We find that the reactive 
NH + NH cross section is driven by a short-ranged collisional mechanism, and its magnitude is only weakly 
dependent on magnetic field strength. Unlike most ultracold reactions observed so far, the NH + NH scattering 
dynamics is non-universal. Our results indicate that chemical reactions can cause more trap loss than spin- 
inelastic NH -I- NH collisions, making molecular evaporative cooling more difficult than previously anticipated. 



The ability to produce and trap molecules at sub-kelvin 
temperatures offers numerous exciting possibilities in chem- 
istry. Recent experiments have demonstrated that ultracold 
chemical reactions can be efficiently manipulated using an 
external electromagnetic field |[l]-[3|], providing new tools to 
control reaction pathways and rate constants. The pronounced 
quantum behavior of ultracold matter may also lead to novel 
phenomena such as "superchemistry", a process in which an 
atomic and molecular Bose-Einstein condensate (BEC) are 
coherently coupled to stimulate chemical reactivity [4]. Up to 
the present, however, studies on ultracold chemical reactions 
have focused only on (bi-)alkali-metal systems with a rather 
limited chemistry. Moreover, most of the observed ultracold 
reactive processes exhibit universal behavior, i.e., the dynam- 
ics are completely determined by long-range interactions [5]. 
Cold reactive collisions in the non-universal regime, as well as 
cold reactions involving non-alkali and open-shell molecules, 
are still largely unexplored. 

Cold chemistry is also relevant in the context of evaporative 
and sympathetic cooling. These second-stage cooling meth- 
ods, which rely on strong elastic collisions between trapped 
particles, represent the final step towards full quantum de- 
generacy and Bose-Einstein condensation |6]. As a rule of 
thumb, the ratio between elastic and non-elastic (inelastic and 
reactive) cross sections should be at least two orders of mag- 
nitude in order for second-stage cooling to succeed. For polar 
molecules, whose properties are expected to find wide appli- 
cations in ultracold physics , cooling into the quantum- 
degenerate regime is yet to be achieved experimentally. Var- 
ious theoretical studies suggest that molecular second-stage 
cooling is feasible |12-19], but these investigations are based 
on the assumption that chemical reactions are strongly sup- 
pressed in a (magnetically trapped) spin-polarized gas ll20ll 
and that the dynamics evolves on a single non-reactive po- 
tential. To our knowledge, this assumption has not yet been 
validated by explicit reactive scattering calculations. 

In this Letter, we present rigorous coupled-channels calcu- 
lations for the reactive NH(X ^S") + NUiX ^S") system in 
the presence of a magnetic field. Ultracold reactive NH + 
NH collisions can yield as many as eight different product 
arrangements |21, 22], making it a versatile system for cold 
chemistry studies. NH has already been cooled to millikelvin 
temperatures and stored in a magnetic trap using Stark decel- 



eration 11231 12411 and buffer gas cooling techniques Ill5ll25l l2^ 
and earlier theoretical studies - based on non-reactive scatter- 
ing calculations - have indicated that NH is a promising can- 



didate for second-stage cooling experiments ||12| - |17LI19LI27[1 . 
Within a magnetic trap, NH is polarized in the low-field- 
seeking state |5nh = 1, -^^Snh — 1)' with Snh denoting the 
total electronic spin and ^/snh projection onto the mag- 
netic field axis. A collision complex of two such molecules 
is in the high-spin quintet \S — 2, Ms — 2) state, with S 
and A/5 referring to the dimer spin quantum numbers. Spin- 
inelastic NH + NH collisions can change either the A/5 quan- 
tum number of the quintet state or the total spin S to produce 
singlet (S — 0) or triplet (S = 1) complexes. Since the S — 
and 1 states are chemically reactive [21i'], S'-changing transi- 
tions may also initiate chemical reaction. 

In order to study the reaction dynamics of cold NH radi- 
cals, we have employed a single-arrangement quantum reac- 
tive scattering method that allows for the calculation of to- 
tal NH + NH reaction probabilities. The cross sections are 
obtained from full quantum scattering calculations on three 
(S — 0, 1, 2) coupled potential energy surfaces. Our results 
indicate that chemical reactions can cause more trap loss than 
inelastic collisions, implying that reactive channels must be 
taken into account when assessing the feasibility of evapora- 
tive and sympathetic cooling. We also find that the total re- 
action probability is strongly dependent on the details of the 
(short-range) interaction potentials, providing one of the first 
examples of non-universal cold chemistry. 

We focus on collisions between two magnetically trapped 
bosonic ^^NH(X '^I]^) molecules and treat the monomers as 
rigid rotors. The coordinate system consists of the intermolec- 
ular vector R, with length R, the polar angles 6a and 9b of the 
monomers, and the dihedral angle (p. The scattering Hamilto- 
nian is given by 
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where /i is the reduced mass of the complex, is the 
angular momentum operator associated with rotation of R, 
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FIG. 1 : Cuts through the quintet (cyan), triplet (magenta), and singlet 
(dark blue) potential energy surfaces of NH(X ^E") - NH(X ^E~) 
for = 180° and 9a ^ 0b = 9. 



Vs{R,OA,OB,(f>) is the potential energy surface for total 
spin S, V;nagn.dip(-R, Sa,Sb) IS the intermolecular magnetic 
dipole interaction between the two monomer triplet spins, and 
Ha and Hb are the Hamiltonians of the individual monomers. 
The latter account for the monomer rotation, intramolecular 
spin-spin coupling, spin-rotation coupling, and Zeeman in- 
teraction [17]. The terms that couple the S = 0,1, and 2 
potentials and that ultimately drive spin-changing processes 
are the intermolecular magnetic dipole interaction and the in- 
tramolecular spin-spin and spin-rotation couplings. 

The NH-NH potential energy surfaces have been obtained 
from high-level ab initio calculations, as described previously 
11281]. Figure [U shows a cut through the three NH-NH poten- 
tials for (p = 180° and 0^ = 0^ = ^- It can be seen that the 
singlet and triplet potentials are strongly attractive at small 
intermolecular distances, which is due to their chemically re- 
active nature. In view of these deep potential energy wells 
we may assume that, once a reactive singlet or triplet NH- 
NH complex is formed, the system readily undergoes exoergic 
chemical rearrangement. For instance, NH + NH may react 
into N2H2 (provided that a third body can dissipate the excess 
kinetic energy) or into a binary product configuration such as 
N2 + H2 02211 . In order to calculate the total reaction probabil- 
ity, we consider only the NH + NH reactant arrangement and 
apply "capture" boundary conditions at short range. That is, 
at a sufficiently small value of the radial coordinate R, we al- 
low flux to disappear into reactive channels. Such an approach 
is commonly used in (reactive) scattering problems involving 
deep potential energy wells ll29ll . We note that not all colli- 
sions on the singlet and triplet potentials are reactive, since 
these surfaces are repulsive for certain geometries. Hence, S- 
changing collisions will not necessarily lead to chemical re- 
action. This is also contained in our boundary conditions (see 
Supplementary Material). ColUsions occurring on the quintet 
potential are entirely non-reactive. 

We have developed a novel reactive scattering algorithm 
based on the renormalized Numerov propagator to calculate 
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FIG. 2: Cross sections for magnetically trapped ^^NH + ^^NH as a 
function of collision energy, calculated for a magnetic field strength 
of 1 G. The solid lines were obtained by including all three potentials 
in the scattering calculations, while the dashed lines were obtained 
using only the non-reactive quintet potential. 



the relevant collision cross sections. Our algorithm amounts 
to solving the coupled-channels equations for two sets of solu- 
tions, and subsequently applying reactive scattering boundary 
conditions to obtain the S'-matrix. Details can be found in the 
Supplementary Material. We used a symmetry-adapted chan- 
nel basis set with even permutation symmetry and even par- 
ity, and a space-fixed total angular momentum projection of 
A4 ~ 2 (see also Supplementary Material). This basis allows 
for s-wave collisions between identical bosonic molecules in 
the same initial quantum state. The radial grid ranged from 
4.5 to 1500 flo, with a minimum of 10 grid points per smallest 
(local) de Broglie wavelength. The cross sections were cal- 
culated for collision energies E of 10^^ to 1 K and magnetic 
field strengths B of lO"! to 10"* G. 

Figure|2]shows the elastic, spin-inelastic, and reactive cross 
sections for two magnetically trapped NH molecules as a 
function of the collision energy for B = 1 G. We find that 
the elastic cross section at small E is constant as a function 
of energy, consistent with Wigner's threshold law for s-wave 
elastic scattering |30], Cciastic E^'"^^^""^^ . Here Ljn and 
Lout denote the partial waves for the incoming and outgoing 
channels, respectively. The inelastic and reactive cross sec- 
tions exhibit E~^l'^ threshold behavior, as expected from the 
^Li„-i/2 i^Qj. exoergic s-wave (Lin = 0) collisions |30i]. 
In order to compare our results with the non-reactive case, we 
have also plotted the cross sections obtained from scattering 
calculations on the non-reactive quintet potential. As can be 
seen in Fig.|2] the inclusion of chemically reactive {S — ^ and 
1) potentials has an almost negligible effect on the elastic and 
spin-inelastic cross sections, confirming our earlier expecta- 
tions reported in Ref. [27]. We thus conclude that most of the 
non-reactive, inelastic trap loss occurs on the quintet surface, 
and that collisions on the S* = and 1 potentials are almost 
100% reactive. 
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FIG. 3: Cross sections for magnetically trapped ^^NH + ^^NH as a 
function of magnetic field strength, calculated at i? = 10~® K. The 
spin-inelastic cross sections obtained with the intermolecular mag- 
netic dipole term switched off, and the reactive cross sections ob- 
tained with the intramolecular spin-spin coupling switched off, are 
also plotted. 
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FIG. 4: Cross sections for magnetically trapped ^^NH -I- ^''NH as 
a function of collision energy, calculated for several magnetic field 
strengths. The elastic cross sections are the same for all magnetic 
fields considered, and the reactive cross sections are the same for 
B — 0,0.1, and 10 G. The results of the universal quantum-defect 
model are also shown. 



Figure 12] also indicates that (for B = 1 G) the elastic -to- 
reactive cross section ratio is much smaller than the elastic-to- 
inelastic ratio, the difference being 2 to 3 orders of magnitude. 
This leads us to reconsider the prospects for molecular evap- 
orative cooling. While previous studies based on non-reactive 
scattering calculations have found that evaporative cooling of 
NH is feasible 1 1/7, 27], our present results indicate that 
chemical reactions can pose a major constraint on the effi- 
ciency of the cooling process. Thus, second-stage cooling of 
NH might be more difficult than previously expected. We will 
return to this topic later, when discussing the effect of uncer- 
tainties in the potentials. 

In Fig. [3] we present the collision cross sections as a func- 
tion of magnetic field for a collision energy of 10~^ K. We 
have already established in previous work lll7[ bill that the 
intermolecular magnetic dipolar term is the dominant source 
of inelastic trap loss at ultralow energies. Indeed, when dis- 
carding the intermolecular magnetic dipole-dipole coupling 
(l^nagn.dip) in our reactive scattering calculations, the spin- 
inelastic cross section decreases by several orders of magni- 
tude. For the reactive cross sections, however, we find that 
the main contribution comes from the intramolecular spin- 
spin coupling (denoted as Vintra-ss in Fig- EJ. The effect of 
the intermolecular term Fmagn.dip on the chemical reactivity is 
almost negligible. This can also be understood by considering 
that the intermolecular spin-spin interaction is long-ranged 
113 in . while a chemical reaction can only proceed when the 
reactants approach each other to a very short distance. Hence, 
the intramolecular spin-spin couplin g, w hich acts through the 
potential anisotropy at short range 33|], plays the most 
important role in enabling the reaction to occur. 

It may also be seen that the reactive cross section at small 
fields (B < 100 G) is constant as a function of B. This 



is essentially a consequence of the £;iin-i/2 fQj- exoer- 
gic collisions: the threshold behavior is determined only by 
the centrifugal barrier in the entrance channel, and the energy 
of the outgoing (reactive) channel has no effect on the cross 
section. Note that a different threshold law applies for short- 
range-induced inelastic processes, for which the centrifugal 
barrier and energy of the outgoing channel do play a role ll34ll . 

Figure |4] compares the scattering results for different mag- 
netic field strengths, including i? = 0, as a function of col- 
lision energy. We find that the reactive cross sections are in- 
dependent of magnetic field for B < 100 G at virtually all 
energies considered, and the {S", E'^'n^i/^} regime extends 
down to i? = 0. The inelastic cross sections, however, are 
different for all B values, and they become constant as a func- 
tion of energy for B = ||36l]. Thus, the reactive and inelastic 
cross sections show fundamentally different behavior in the 
ultracold regime. 

Let us now consider ultracold NH + NH scattering in the 
universal limit. Knowledge of the degree of universality is 
useful to understand and interpret the (generally complex) 
full-dimensional reaction dynamics in terms of simple few- 
parameter models. In order to establish whether the ultracold 
chemistry of NH + NH is universal, we compare our scatter- 
ing results with the single-channel quantum-defect model of 
Idziaszek and Julienne |5]. Details of this model are given in 
the Supplementary Material. In the universal limit, all scat- 
tering flux that reaches the short range disappears into reac- 
tive channels. Figure |4] compares the results of the universal 
quantum-defect model with those of the numerical coupled- 
channels calculations. In the s-wave regime (E < 10~^ K), 
the elastic cross sections are underestimated by the universal 
model, while the reactive cross sections are clearly overesti- 
mated. The ratio between the universal elastic and reactive 
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FIG. 5: Cross sections for magnetically trapped ^^NH + ^''NH as 
a function of the scaling factor A, calculated for two magnetic field 
strengths (B = 1 and 100 G) at E = W'^ K. The results of the 
universal quantum-defect model are also shown. 



cross sections differs by more than one order of magnitude 
(a factor of w 26) from the numerical data. These differ- 
ences suggest that the NH + NH reaction dynamics is non- 
universal, because a significant fraction of the incident flux 
is reflected at short range. It should be noted, however, that 
the quantum-defect results apply to single-channel scattering 
on a single isotropic potential, while the numerical data have 
been obtained from multi-channel calculations on three cou- 
pled anisotropic potential energy surfaces. Nevertheless, the 
collision dynamics in the ultracold regime is s-wave domi- 
nated, and the effective NH-NH long-range potential - which 
is the same for all three spin states - is governed mainly by 
the isotropic interaction. Thus, the single-channel quantum- 
defect model should be applicable to NH-NH in the universal 
regime. 

An alternative way to establish the degree of universality 
is to test the effect of small modifications in the (short-range) 
potentials. If the scattering is universal, the cross sections are 
completely determined by the long-range features of the in- 
teraction potentials. For instance, a scaling of the potentials 
by a factor of A (or, equivalently, a scaling of the reduced 
mass by A) should change the universal elastic cross section 
by A^/^. Figure |5] shows the universal quantum-defect results 
as a function of the scaling parameter A (0.9 < A < 1.1) for 
E = 10~^ K. The corresponding numerical cross sections, 
obtained by reduced-mass scaling, are also shown for B = 1 
and 100 G. It is evident that the numerical results are highly 
sensitive to the details of the potentials, and that the univer- 
sal model is inaccurate for all values of A. In fact, the reso- 
nance features in the numerical cross sections are signatures 
of non-universal behavior [SS'l and highlight the importance of 
short-range physics in the dynamics. We thus conclude that 
the scattering properties of magnetically trapped NH cannot 
be captured in a universal model. 

As discussed extensively in Refs. lIB. 17, 27], the A-scaling 



approach also provides a means to sample the effects of uncer- 
tainties in the interaction potentials. Although the potentials 
have been obtained from state-of-the-art methods, there is a 
small remaining inaccuracy in the ab initio data which gives 
rise to an uncertainty in the calculated cross sections. This 
also carries implications for the prospects for molecular evap- 
orative cooling. More specifically, one should evaluate the 
elastic-to-reactive cross section ratio for all relevant A values 
to obtain a realistic estimate of the cooling efficiency. In the 
case of A = 1, we find that reactive NH + NH collisions are 
more probable than elastic ones (cf. Figs. |3] and |4]i, while for 
A « 0.95 and A ~ 1.03 the elastic cross sections are about one 
order of magnitude larger than the reactive ones. These results 
suggest that evaporative cooling of magnetically trapped NH 
might still be feasible, but the probability of success is signif- 
icantly smaller than estimated earlier from non-reactive scat- 
tering calculations. We note that these findings are still valid 
when the size of the channel basis set is increased, as detailed 
in the Supplementary Material. 

It can also be seen in Fig.|5]that the inelastic cross sections 
change rather dramatically from S = 1 to 100 G, while the 
reactive cross sections show only a weak dependence on mag- 
netic field. Nevertheless, for certain values of A, the reactivity 
can increase by almost one order of magnitude as the mag- 
netic field strength is changed. Thus, it may be possible to 
control the NH + NH reaction rate by means of an external 
field. For most A values, however, the reactive cross sections 
show 5° behavior and magnetic field control will not be pos- 
sible. We note that the final product-state distribution might 
be more sensitive to the magnetic field strength than the to- 
tal reaction probability, but our scattering method does not 
allow the calculation of product-state-resolved reaction cross 
sections. 

In summary, we have presented the first rigorous quan- 
tum scattering study of ultracold reactive molecular colli- 
sions in the presence of a magnetic field. Our results illus- 
trate the importance of chemical reactions in a magnetically 
trapped molecular gas, and call to reconsider the prospects for 
molecular evaporative cooling and magnetically-controlled 
cold chemistry. This work may serve as a benchmark for other 
ultracold paramagnetic and dipolar molecules, the dynamics 
of which are still virtually unexplored. 
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